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X-ray variability in a deep, flux limited sample of QSOs 



o 
o 
o 

(N 

Oh- 
< 

(N 
> 
00 

O 

(N 

ON 

On 



9 L,: 
6 



O. Almaini, 1 A. Lawrence, 1 T. Shanks, 2 A. Edge, 2 B.J. Boyle, 3 I. Georgantopoulos, 4 
K.F. Gunn, 5 G.C. Stewart, 6 and R.E. Griffiths 7 

1 Institute for Astronomy, Royal Observatory, University of Edinburgh, Blackford Hill, Edinburgh EH9 3HJ 

2 Department of Physics, University of Durham, South Road, Durham DH1 3LE 

3 Anglo-Australian Observatory, PO Box 296, Epping NSW 2121, Australia 

4 National Observatory of Athens, I. Metaxa & B. Pavlou, Palaia Penteli, 15236, Athens, Greece 

5 Department of Physics & Astronomy, University of Southampton, University Road, Southampton SO 17 IB] 
B Department of Physics and Astronomy, The University of Leicester, Leicester LEI 7RH 

7 Department of Physics, Carnegie Mellon University, Wean Hall, 5000 Forbes Ave., Pittsburgh, PA 15213, USA 



MNRAS in press 



ABSTRACT 

We present an analysis of X-ray variability in a flux limited sample of QSOs. Selected from 
our deep ROSAT survey, these QSOs span a wide range in redshift (0.1 < z < 3.2) and are 
typically very faint, so we have developed a method to constrain the amplitude of variability 
in ensembles of low signal-to-noise light curves. We find evidence for trends in this variability 
amplitude with both redshift and luminosity. The mean variability amplitude declines sharply 
with luminosity, as seen in local AGN, but with some suggestion of an upturn for the most 
powerful sources. We find tentative evidence that this is caused by redshift evolution, since the 
high redshift QSOs (z > 0.5) do not show the anti-correlation with luminosity seen in local 
AGN. We speculate on the implications of these results for physical models of AGN and their 
evolution. Finally, we find evidence for X-ray variability in an object classified as a narrow 
emission-line galaxy, suggesting the presence of an AGN. 
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1 INTRODUCTION 

The X-ray flux from AGN exhibits variability on shorter timescales 
than any other waveband, indicating that this emission occurs in 
the innermost region of the central engine. Detailed studies of X- 
ray variability can therefore act as important probes of the AGN 
phenomenon, potentially yielding information on size scales unre- 
solvable by any other means. 

While previous observations had revealed that X-ray variabil- 
ity was very common in AGN (see Mushotzky, Done & Pounds 
1993 for a review), the most significant advances were made with 
the 'long look' observations taken with the EXOSAT observa- 
tory. The highly eccentric orbit of this satellite combined with its 
much improved sensitivity allowed continuous observations of a 
small number of local AGN with high signal to noise. The light 
curves appeared very random in nature. Fourier and spectral analy- 
ses revealed no preferred timescales or periodicities. Instead these 
sources showed power smoothly distributed over a wide range of 
frequencies, generally with power spectra of the form P(f) oc 
f~ 1 - 5 (Lawrence & Papadakis 1993, Green et al 1993). Such light 
curves are described as red noise, i.e. scale free but with more 
power at low frequencies. The power spectra are eventually ex- 
pected to flatten from a power-law form at the lowest frequencies 
below ~ 10~ B Hz (Edelson et al 1999). Similarly, at high frequen- 



cies we expect a cutoff corresponding to the smallest scales in the 
emitting region (e.g. the inner region of the accretion disk) but so 
far no such steepening has been observed. 

Interesting results were obtained by comparing the variabil- 
ity in different AGN, suggesting that more luminous sources vary 
with lower amplitude. Barr & Mushotzky (1986) studied the short- 
est timescale required for a source to double in flux and observed an 
increase in the doubling timescale with luminosity. Lawrence & Pa- 
padakis (1993) and Green et al (1993) studied the power spectra of 
the EXOSAT AGN and found that the power at a given frequency 



decreases with luminosity with the form P oc L, 



-(0.6±0.1) 



This 



result has since been confirmed by observations with the ASCA 
satellite, albeit with overlapping samples of AGN (Nandra et al 
1997, Turner et al 1999). Such trends might be explained if more 
luminous sources are physically larger in size. However, these stud- 
ies were conducted on a small number of bright AGN at low red- 
shift (z < 0.1). Since QSO activity is known to have peaked at 
much earlier epochs (z ~ 2) it is now vital to establish whether 
the trends observed in these local AGN apply to the QSO popu- 
lation as a whole and to test for evolutionary effects. In this work 
we present an analysis of a flux limited sample of 86 QSOs de- 
tected as part of the Deep ROSAT Survey (Georgantopoulos et al 
1996). Spanning a wide range in redshift (0.1 < z < 3.2), the 
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same sample has been used to define the form and evolution of the 
QSO luminosity function (Boyle et al 1994). Most of the sources 
are detections at the faintest possible X-ray flux and hence the light 
curves are generally of low signal to noise. Nevertheless, we have 
developed methods for constraining the level of intrinsic variability 
in each source. By combining the results from many QSOs we can 
search for trends in the typical variability amplitude with redshift 
and luminosity. Finally we study the variability of X-ray emitting 
galaxies to test the hypothesis that these contain obscured AGN. 



2 THE SAMPLE 

Our sample consists of 86 QSOs identified during the Deep ROSAT 
Survey of Shanks et al (in preparation), which reached an X-ray 
flux limit of ~ 4 x l(n 15 erg s^cm -2 (0.5-2.0 keV). A descrip- 
tion of the X-ray source detection and optical identification is given 
in Georgantopoulos et al (1996). Details of the QSO redshift dis- 
tribution and the X-ray luminosity function and its evolution are 
given in Boyle et al (1994). 

The survey uses 7 deep ROSAT PSPC pointings, with integra- 
tion times of ~ 30 — 80ks per field, spread over periods ranging 
from 2-14 days. Periodic gaps in the light curves are caused by 
the orbit of the satellite, with larger gaps occurring during phases 
of high instrument background, such as when the satellite passes 
through the South Atlantic Anomaly. During each 96 minute orbit, 
a source is typically visible for 1000 — 2000 seconds. Given the 
faint nature of the X-ray sources, light curves were first binned on 
these periodic orbits, giving 16 — 26 observation slots (depending 
on the field). The fainter sources are then re-binned to allow mean- 
ingful Gaussian statistics. Using the mean intensity of the source, 
if the smallest bin is expected to give fewer than 20 photons it is 
merged with its nearest neighbour. This process is repeated until 
I At > 20 photons per bin, leaving 86 QSOs with an absolute min- 
imum of at least 2 time bins. We note that this re-binning process 
will inevitably smear out high frequency variability, and discuss 
this effect further in Section 3 and Appendix A. 

Basic data reduction was carried out using the Asterix X- 
ray data processing package. For each source, depending on the 
off-axis radius, a source box radius was chosen to enclose 90 per 
cent of the X-ray photons. Assuming a mean QSO spectral index of 
r = 2.2 (Almaini et al 1996), the radius of this circle varies from 
24.6" on axis to 56.6" at a maximum radius of 18 arcmin. Data 
from periods of high particle background are excluded from the 
analysis. This removes ~ 10 per cent of the data when the Mas- 
ter Veto Rate rises above 170 counts s _1 (Plucinsky et al 1993). 
Considerable care was taken in choosing areas for background sub- 
traction. By good fortune, none of the fields showed evidence for 
significant background gradients (generally caused by irregulari- 
ties in the galactic background or contamination from solar scat- 
tered X-rays). Typically 10 circular areas of 4-6 arcmin radius are 
then chosen across each field from source free regions to perform 
the background subtraction. The background box closest to each 
source was used, with minor vignetting corrections (< 5 per cent) 
to correct for the differing off axis radii of source and background 
boxes. 

Problems may also conceivably arise if the background shows 
differing variability across the field. Although care is taken to use 
background regions close to each source, differential variability 
may nevertheless introduce extra variance into the light curves 
above the Poisson noise and bias the results. Some of the fields do 
indeed show evidence for variability at a low level (51/ 1 < 10 per 



Figure 1. Light curve and likelihood function (of arbitrary normalisation) 
for the QSO GSGP4X:015 
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cent), probably due to contamination from solar scattered X-rays, 
but this is consistent across the detector. Such minor effects will 
therefore be removed in background subtraction. Finally, to check 
that any variability was real and not due to a flickering response 
in the PSPC, we compared the light curves of brighter QSOs with 
each other and with background regions. No significant correlations 
could be seen, suggesting that the bulk of the variability is indeed 
real and not caused by mysterious instrumental effects. 



3 MEASURING THE AMPLITUDE OF QSO 
VARIABILITY 

We now develop a method for measuring the amplitude of intrinsic 
fluctuations in light curves of low signal to noise. Although many 
will be too faint to enable a significant detection of variability, we 
can nevertheless place upper limits on the amplitude of fluctuations 
which could be present. In order to combine information from faint 
light curves, we will also extend this method to measure the typical 
variability in an ensemble of QSOs (e.g. within a given luminosity 
or redshift range). 

Unlike the well studied EXOSAT light curves, the QSOs 
studied here are too faint to obtain variability amplitudes by de- 
tailed power spectrum analysis. A simpler method is to estimate 
the intrinsic variance in the light curve. Since this quantity is pro- 
portional to the amplitude of the power spectrum (see Appendix 
A) it will also allow comparisons with previous work. Nandra 
et al (1997) used the root mean square variation in their ASCA 
lightcurves, subtracting off a contribution due to noise. An intrin- 
sic assumption in this method is that each point in the light curve 
has equal weight, which will break down in cases where the errors 
differ significantly from point to point. We therefore develop a new 
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Figure 2. Light curve and likelihood function (of arbitrary normalisation) 
for the QSO GSGP4X:040 



using Bayes' Theorem to obtain the probability distribution for oq 
given our measurements: 
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p(a Q \xi,ai) = p(x i \a i ,a Q ) 
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where L (<jq \xi, Oi) is the likelihood function for the parame- 
ter oq given the data. This general form for the likelihood function 
can be calculated if one assumes a Bayesian prior distribution for 
oq and Xi. In the simplest case of a uniform prior one obtains: 



L(<j Q \xi,(Ji) oc p(xi\ai,aQ) 
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By differentiating the above, the maximum likelihood esti- 
mate for <tq can be shown to satisfy the following, which (for a 
uniform prior) is mathematically identical to a least \ 2 solution: 



» {( Xl -x)*-(ai + a 2 Q )} 

4f K 2 + ^) 2 



(7) 



In the case of zero measurement errors, this reduces to the 
standard form: 



maximum likelihood technique to extract the best estimate of the 
intrinsic (noise subtracted) variance given a set of data values Xi 
with differing measurement errors Oi. 

In order to compare objects of different flux in a meaningful 
way, all light curves are first divided by their mean flux. Hence we 
are measuring the amplitude of fractional variability. When an en- 
semble of QSOs is considered (e.g. all QSOs within a given redshift 
range) we treat each light curve as a 'snapshot' from the population 
and assume that the ensemble all have the same underlying distri- 
bution in 7/7. By definition the variance does not depend on the 
ordering of the points and hence we can combine the light curves 
and measure the variance in the resulting distribution. We now de- 
scribe a method for separating the intrinsic variance from the scatter 
due to measurement errors. 



3.1 A maximum likelihood estimator 

The variance in a light curve has two components: first there will 
be fluctuations due to noise (cr 2 oise ) and secondly there may be in- 
trinsic variations from the QSO or ensemble of QSOs (<Tq). Hence: 
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If we assume Gaussian statistics, for a light curve with mean 
x and measurement errors ou the probability density for obtaining 
the N data values Xi is given by: 



p(Xi\<Ti,<T Q ) = Y\ 
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This is simply a product of N Gaussian functions representing 
the probability distribution for each bin. We may turn this around 



lim <jq — y~^(a:j — x) 2 /N 



(8) 



In the case of identical measurement errors (a m = constant) 
this reduces to the form used by Nandra et al (1997): 



lim <tq = { (xi - x) 2 - o 2 m ) /N 



(9) 



In general however, Equation 7 can only be solved analyti- 
cally for the case where N < 2. For a larger number of bins we 
must find the maximum likelihood estimate numerically. Two con- 
trasting examples of light curves and their likelihood functions are 
given in Figures 1 and 2, from QSOs of high and low flux respec- 
tively. From these, maximum likelihood estimates for ctq can be 
evaluated with suitable error bounds. The maximum likelihood es- 
timate for ctq is obtained by locating the peak of the likelihood 
curve, while the error bounds are obtained in the standard way by 
assuming that the likelihood curve has a Gaussian shape. 

This method can either be applied to individual light curves or 
an ensemble. In comparing values of oq, however, it is important 
to be aware of the effects of different total integration times and 
sampling. In addition, one might expect cosmological time dilation 
to modify the results. A detailed examination of these effects can be 
found in Appendices A & B. We find that one can apply a correc- 
tion to allow for these changes, but only under the assumption of a 
particular underlying power spectrum. We therefore carried out our 
analysis with and without these timescale corrections, but overall 
we found that they have no significant effect on the results. 



4 THE RESULTS 
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Figure 3. Light curves of the 12 QSOs showing evidence for variability 
with > 95 per cent significance. 



SGP3X:035 



Figure 4. Showing the maximum likelihood estimates for the intrinsic vari- 
ability as a function of flux for the 86 QSOs. Timescale corrections have 
not been applied. The point size is proportional to the object flux, to aid in 
distinguishing objects of higher S/N in Figures 5, 6, 8 and 9. 




0.1-2keV flux, erg s~ l crrr 2 



4.1 Detecting individually variable QSOs 

Before applying the maximum likelihood technique, a simple x 2 
test was performed on each light curve to test for variability against 
the null hypothesis that the flux remains constant: 



X 



\- [In ~ {In}] 2 



(10) 
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The results show that only 12 of the 86 QSOs show individual 
evidence for variability with a significance > 95 per cent, of which 
4 are significant at > 99 per cent. The light curves for these QSOs 
are displayed in Figure 3. As expected, the fraction of significantly 
variable objects increases at higher flux. Of the 10 brightest QSOs 
for example, 5 show variability at > 95 per cent significance, with 
a further 2 showing variability at a significance of 90 per cent, sug- 
gesting that a large fraction of all QSOs would show variability 
given sufficient signal-to-noise. 

It is important to emphasise that the non-detection of variabil- 
ity in the remaining QSOs does not imply that no variability is oc- 
curring. Any intrinsic variations which are present are simply over- 
whelmed by photon noise. The detection of significant variability 
depends on the amplitude of intrinsic fluctuations and on the signal 
to noise. Since we can measure the photon noise, it should there- 
fore be possible to at least place an upper limit on the amplitude of 
intrinsic fluctuations in every QSO. Only then can we make quan- 
titative statements about the "typical" variability in the population. 



Time (seconds) 



4.2 The maximum likelihood results 

Numerically solving Equation 7 for each of the 86 QSOs, the maxi- 
mum likelihood estimates for the intrinsic variability amplitude are 
displayed in Figure 4 (with 68 per cent confidence regions), plot- 
ting a as a function of flux. The timescale corrections outlined in 
Appendix A are not applied to these raw data points, but these have 
no major effect on the results (see below). As expected, for most 
faint QSOs individually we can only estimate upper limits for the 
variability amplitude, but nevertheless the brighter QSOs and the 



© 0000 RAS, MNRAS 000, 000-000 



X-ray variability in a deep, flux limited sample of QSOs 5 



Figure 5. (a) Maximum likelihood estimates for the variability amplitude 
as a function of luminosity for the 86 QSOs. The point size is proportional 
to the object flux (see Figure 4) in order to highlight objects with higher 
S/N. In (b) we display the results in ensemble form, with unfilled squares 
showing the minor effect of not applying timescale corrections. 



Figure 6. (a) Maximum likelihood estimates for the variability amplitude as 
a function of redshift for the 86 QSOs. The point size is proportional to the 
object flux (see Figure 4) in order to highlight objects with higher S/N. In 
(b) we display the results in ensemble form, with unfilled squares showing 
the minor effect of not applying timescale corrections. 





fainter objects combined would seem to suggest a typical variabil- 
ity amplitude significantly above zero (with a ~ 0.2). 

In Figures 5(a) & 6(a) we display the variability amplitude as a 
function of luminosity and redshift. We split the data into 6 equally 
spaced luminosity and redshift bins and measure the ensemble vari- 
ance, applying the same maximum likelihood technique to the com- 
bined light curves (Figure 5(b) and 6(b)). The ensemble results are 
displayed with and without the timescale corrections, which are de- 
pendent on the assumed power spectrum. Although the net effect of 
these corrections is to boost a upwards, they clearly have no sig- 
nificant effect on the results. 

We find a clear detection of variability in the ensemble of 
QSOs, with evidence for trends with luminosity and redshift. A 
X 2 test was applied to test the null hypothesis that the scatter in 
the 5 points on Figures 5(b) and 6(b) is consistent with their error 



bounds. For both the redshift and luminosity data, we can reject the 
null hypothesis with > 99 per cent confidence. 

4.3 Correlations With Redshift And Luminosity 

Previous studies of very local (z < 0.1) AGN found significant ev- 
idence for a relationship variability amplitude and luminosity, pa- 
rameterised by the form a 2 oc Z /a r ( ' ' 6±0 ' 1 - ) (Lawrence & Papadakis 
1993, Green & McHardy 1993, Nandra et al 1997). In Figure 7 we 
compare our ensemble QSO results with the work of Papadakis 
& Lawrence (1993) and Nandra et al (1997) (although we note 
that there is considerable overlap between these local samples). 
Similar results have recently been found by Turner et al (1999). 
To allow a meaningful comparison we convert luminosities to our 
0.1 — 2 keV band (assuming a spectral shape of the form a = 1.0 
above 1 keV and a = 1.5 below 1 keV). Values of a were cor- 
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Figure 7. A comparison of the ROSAT 'results with the findings of Lawrence 
& Papadakis (1993) (circles) and Nandra et al (1997) (triangles). The best 
fitting power-law of the form a 2 oc L x 0,(i is shown (not fitted to the 
ROSAT QSOs). To allow a crude comparison, values of a have been nor- 
malised to timescales below 1 week using a temporal power spectrum of the 
form / ~~ 1 ' 5 and luminosities have been converted to the . 1 — 2 keV band. 

rected to timescales of 1 week (assuming a power spectrum of 
the form P(f) oc / _1 ). Given the uncertainties in these conver- 
sions, and the fact that previous variability studies were conducted 
at higher X-ray energies, we will be cautious in making direct com- 
parisons. However, we note that the amplitudes found in our QSOs 
are broadly consistent with these local AGN. This would support 
the hypothesis that the driving force in any variability trend is ac- 
tually the QSO luminosity, rather than the redshift. Formally, how- 
ever, a power-law gives a very poor fit to the data (see Table 1). 

An intriguing result is the apparent minimum at L x ~ 5 x 
10 44 erg s _1 , with an upturn towards higher luminosities. No evi- 
dence for such an upturn has been seen among local AGN. To allow 
a fairer comparison, we split our sample into 'local' and 'distant' 
QSOs, dividing at z = 0.5 (Figures 8,9). We find that while the lo- 
cal QSOs show a clear anti-correlation with luminosity, consistent 
with previous studies, the QSOs at higher redshift show no evi- 
dence for this trend whatsoever. To parameterise this difference, we 
perform power-law fits to these distributions (see Table 1 ; similar 
results are obtained if we fit to the binned, ensemble data). Unlike 
the full sample, we find that a power-law can describe these sep- 
arate distributions well (Xred ~ 1) an d we obtain good fits with 
a 2 oc " ,±0 ' 2 ' for the low redshift QSOs but a very different 
form for the higher redshift objects, with a 2 oc £+(°- 6±0 - 3 ) (j e 
a possible positive correlation with luminosity). There is a strong 
degeneracy between luminosity and redshift, however, and it would 
clearly be desirable to obtain a larger sampling of the L — z plane. 
If one restricts the comparison to only the 3 overlapping luminos- 
ity bins, for example, a x 2 test rejects the hypothesis that the two 
distributions are the same at only the 95 per cent confidence level. 
We therefore urge caution in over-interpreting these results, but we 
will explore some of the possible consequences in Section 5.4. 



Table 1. Showing power-law fits to the distribution of a 2 vs luminosity for 
the ROSAT QSOs. 

Data Best fit X 2 red 

All QSOs a 2 ocL- (0 ' 8±0 ' 2) 2.5 
Z <0.5 ^£-(1.5*0.2) 12 

z > 0.5 a 2 oc L+(°- 6±0 - 3 ) 0.8 



5 CONSTRAINING PHYSICAL MODELS 

Although we find some evidence for an upturn in the a/L x relation 
at high luminosities, it is now well established that a strong anti- 
correlation exists between these two parameters, at least for local 
AGN. We will now explore some possible physical implications of 
this trend, returning to the possible consequences of an upturn in 
Section 5.4. 

The most plausible explanation for the observed anti- 
correlation is that high luminosity sources have larger emitting re- 
gions. It is currently unclear whether this is due to a single, coher- 
ent but erratically varying source or a large number of independent 
flaring regions (Abramowicz et al 1991). The incoherence in the 
light curves of EX OS AT AGN would suggest a stochastic random 
process governed by a large number of variables (Krolik, Done 

6 Madejski 1993), as would the lack of any obvious timescale 
(Lawrence et al 1987). We now investigate how the observed re- 
lation a 2 oc L x ( fi±01 ) can be used to constrain two classes of 
models for AGN variability and in turn discuss the implications for 
QSO evolutionary models. 

5.1 The single coherent oscillator 

In the case where the X-ray variability is caused by global, coher- 
ent changes in the emitting region, it is reasonable to assume that 
the luminosity of the source will be related to its physical size. The 
geometry of the emitting region is unclear however, and could be 
spherical, disk shaped etc., and hence we assume a general relation- 
ship between the X-ray luminosity and the characteristic radius, R: 

L x xR (11) 

Over the range of frequencies we are interested in, we will as- 
sume a temporal power spectrum of the form P(f) oc f~ a (where 
a ~ 1.5) as observed in local AGN (Lawrence & Papadakis 1993). 
Assuming self similar scaling, timescales will be scaled in direct 
proportion to the size of the source, and hence in larger sources the 
temporal power spectrum will simply be shifted in the frequency 
direction. Self similar scaling requires conservation of the total in- 
tegrated fractional power, but the power measured over a narrow 
frequency range will change. By integrating the new power spec- 
trum over the observed frequency range one can relate the observed 
variance to the size of the source: 



= / P(f)df (12) 
Jfi 

oc fl~" (/i«/ 2 ,q/1) (13) 

oc R 1 -" (14) 

oc (15) 



Thus one can readily reproduce the observed anti-correlation 
with luminosity (a 2 oc L~ {0 - 6±01) ) given the measured power 
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Figure 8. (a) Variability amplitude as a function of luminosity for the low Figure 9. (a) Variability amplitude as a function of luminosity for the high 

redshift (z < 0.5) QSOs. In (b) we display the results in ensemble form. redshift (z > 0.5) QSOs. In (b) we display the results in ensemble form. 





spectrum slope (a ~ 1.5) if the geometry of the source is such 
that L x oc ii ' 7-1 ' . This would be consistent with a ring of X-ray 
emitting material rather than a disk or spherical cloud. We note, 
however, that alternative geometries could also reproduce the ob- 
servations if we allow a radial decline in the X-ray emissivity (e.g. 
free-free emission in a spherical cloud of material if the electron 
density falls off at large radii). In conclusion therefore, the observed 
anti-correlation can be reproduced in this model and in principle 
could be used to constrain the geometry of the emitting material. 



5.2 The independent hot spot model 

In this class of model, the variability is caused by independent flar- 
ing events. More luminous sources will simply have more flaring 
regions and hence the r.m.s. fluctuations will be smaller. If the 
time-averaged luminosity from each flare is I, the typical number 
of flares at a given time will be related to the total X-ray luminosity, 
L x by: 



N = L x /l (16) 

Thus the fractional r.m.s. variation on the number of flares at 
a given time will be: 

a = VN/N oc l-^J (17) 

If the flares are the same size in all QSOs, this would lead to 
a 2 oc LjT 1 , which is somewhat steeper than the observed correla- 
tion. This could be overcome if higher luminosity objects had more 
luminous flares, such that I oc L x ' 5 . Hence this model can also 
explain the observed anti-correlation. 

5.3 Constraining models of QSO evolution 

As outlined in Boyle et al (1987) there are at least two possible 
models which can explain the strong evolution in the QSO lumi- 
nosity function: 

(i) The evolution of successive generations of short lived (~ 
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Figure 10. Comparison with the predictions of a long-lived single oscillator 
model. The local trend between variability amplitude and source luminosity 
(a 2 oc Lx 0,6 ; solid line - see Figure 8) is evolved to predict the relation- 
ship at high redshift (dashed lines). For comparison, the ROSAT QSOs are 
superimposed, split into 3 redshifts ranges. 

10 8 yr) QSOs in which the mean intensity of QSO cycles dimin- 
ishes. In this scenario most galaxies have been through a quasar 
epoch at least once. 

(ii) The evolution of individual, long-lived QSOs which uniformly 
dim over time. In this scenario only ~ 1% of galaxies have ever 
undergone quasar activity, but would develop very massive black 
holes. 

In the first model, we expect a typical QSO to harbour a black 
hole of mass ~ 1 x 10 8 M Q regardless of redshift. Black holes 
of this size allow the theoretically favoured accretion at Eddington 
limit (Rees 1984) for a QSO at peak activity. Thus the size of a 
typical emitting region is expected to be very similar at all redshifts. 

In the more controversial long-lived model, QSOs have poured 
out a vast amount of energy over most of the age of the universe 
and should therefore develop extremely massive black holes by the 
associated accretion. This can be calculated by integrating the total 
energy emitted by the QSO from it's formation until the present 
day (full details given in Appendix C). If we assume an accretion 
efficiency e = 0.1 after Rees (1984), an initial mass given by the 
Eddington limit and a starting redshift of z = 3 then we predict a 
factor ~ increase in black hole mass from z — 3 to z = 2, 

with a further growth by a factor of ~ 3/igQ 1 from z = 2 to the 
present day (Figure CI). 

In principle therefore, since variability traces the size scales 
in the central regions of an AGN, one might be able to distinguish 
between these models. The consequences depend on the variability 
model: 

Single Coherent Oscillator: In this case, the amplitude of 
variability is directly related to the physical size of the source 
(Equation 14). The X-ray emission from AGN is generally thought 
to arise a few Schwarzschild radii from the central black hole, and 



in this metric the central mass is directly proportional to the size 
scale. Hence one can relate the amplitude of variability to the cen- 
tral black hole mass, Mbh'- 

a 2 oc R 1 ' oc M B £ 5 (18) 

In the long-lived model this would lead to dramatic changes 
in the variability properties of QSOs with redshift (at a fixed lu- 
minosity). First, as the black hole mass of an individual QSO 
decreases towards higher redshift, the variability will increase as 
a oc Mbh ' 25 (most dramatically at z > 2). In addition, because 
of the evolution in the luminosity function (in this model due to 
the fading of individual QSOs) one should compare objects which 
have followed a similar evolutionary path (i.e. the same starting lu- 
minosity). The observed evolution in the QSO luminosity function 
(Boyle et al. 1994) would imply a decline in luminosity by a factor 
of ~ 30 from z = 2 until the present day. Combining both of these 
effects, the predicted changes are striking (see Figure 10). 

The model predictions are, of course, strongly dependent on 
the assumed accretion efficiencies, redshift of formation etc., but 
the generic prediction is a rise in the mean variability (at fixed lumi- 
nosity) towards higher redshift. Our data do not appear grossly in- 
consistent with these trends, particularly at the highest luminosities 
where the high redshift QSOs do indeed show evidence for an up- 
turn (see also Figures 8 & 9). At present, therefore, we cannot, con- 
clusively rule out a long-lived QSO model on the basis of variability 
evolution. Further measurements are clearly required. In particular, 
examining variability properties at high redshifts (z ~ 3) would 
provide a conclusive test, where the differences are most striking. 
This will be possible with Chandra and XMM. 

In the alternative model, most galaxies have been through a 
quasar phase (for perhaps 10 8 years) and we expect a similar range 
of black hole masses at any given epoch. This model is also con- 
sistent with the current data, although we note that this would not 
explain the apparent upturn in the a/L x relationship for higher red- 
shift AGN. 

Multiple hotspot model: In this model, the variability ampli- 
tude is determined by the number of independent flares. Hence pre- 
dictions relating to evolutionary models will depend on the physics 
of the flaring mechanism. If the number of flares is strictly governed 
by only the source luminosity, then we expect no difference be- 
tween the long and short lived evolution models. However is seems 
plausible that a physically larger emitting region will give rise to 
a greater number of disconnected, independent emitting regions. 
Hence the similarity in the variability trends at high and low red- 
shift may also be sufficient to rule out a long-lived evolutionary 
model. 

5.4 Explaining the high luminosity upturn 

In Section 4 we found evidence that the local anti-correlation be- 
tween variability and luminosity does not hold at the highest lumi- 
nosities. Although further observations are clearly required to test 
the reality of this feature, we will nevertheless speculate on some 
possible explanations. 

Our results tentatively suggest that this feature is caused by 
redshift evolution, since while our low redshift QSOs show the 
same anti-correlation between a and L x as previous studies, the 
higher redshift sources show very different behaviour. A simple in- 
terpretation is that earlier QSOs (at a given luminosity) were pow- 
ered by less massive black holes compared to local AGN (or alter- 
natively were accreting more efficiently). 

On the other hand, if luminosity is the driving parameter 
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Figure 11. Light curve of the variable NELG GSGP4X:091 

(rather than redshift evolution) it could reflect an upper limit to 
the maximum black hole size. Assuming accretion at close to the 
Eddington limit, the upturn X-ray luminosity (5 x 10 44 erg s~ 4 ) 
corresponds to a mass of ~ 5 x 10 7 M Q . This is strikingly sim- 
ilar to the mass at which the space density of active black holes 
turns over sharply (e.g. see Padovani, Burg & Edelson 1990). This 
could, in principle, lead to a levelling off in the a vs. L x relation- 
ship if we assume (not unreasonably) that there is some scatter in 
the luminosity for a given black hole size. 



6 VARIABILITY OF X-RAY LUMINOUS GALAXIES 

In recent years, is has become clear from deep X-ray surveys 
that a population of X-ray luminous narrow emission-line galaxies 
(NELGs) might provide an explanation for the origin of the cosmic 
X-ray background (Boyle et al. 1995, Roche et al. 1995, Carballo et 
al. 1995, Almaini et al 1997, McHardy et al 1998). These galaxies 
are typically 100 times brighter than normal field galaxies, and ap- 
pear to have harder X-ray spectra than QSOs (Almaini et al 1996). 
The nature of this X-ray emission is still controversial, but the pres- 
ence of massive starburst activity and/or obscured AGN has been 
postulated. 

Here we investigate the X-ray variability of a small number 
of these galaxies from our Deep ROSAT Survey. Unfortunately 
most of these galaxies are too faint to enable a meaningful variabil- 
ity analysis. Only 6 of the 23 X-ray luminous galaxies presented 
in Almaini et al (1996) meet the selection criteria outlined in Sec- 
tion 2. A x 2 test was performed on these light curves. The null 
hypothesis of constant flux can only be rejected in the galaxy with 
highest flux, namely GSGP4X:091, with a significance of > 95 
per cent. The light curve is displayed in Figure 1 1 . Such variability 
clearly favours an AGN origin for the X-ray emission. Applying the 
maximum likelihood technique, we obtain variability amplitudes 
a — 0.31±o;i|, consistent with the typical fluctuations seen in low 
luminosity AGN (see Figure 5). Combining the light curves of the 
5 fainter galaxies, we do not detect a variability amplitude signifi- 
cantly above zero. Formally we obtain a 68 per cent upper limit of 
a < 0.19, from which we certainly cannot rule out the presence of 
AGN-like variability. 



7 CONCLUSIONS 

We have developed a technique for constraining the amplitude of 
variability in light curves of low signal to noise. By treating each 



(normalised) light curve as a 'snapshot' from an ensemble, one can 
obtain a maximum likelihood estimate for the typical variability in 
the population. 

We apply this technique to the light curves of 86 QSOs ob- 
tained as part of the Deep ROSAT Survey. This represents the 
first large study of X-ray variability in QSOs over a range in red- 
shift. On timescales of several days, the amplitudes of variability 
are very similar to those found in very local AGN. We find evi- 
dence for trends in the ensemble variability amplitude with both 
redshift and luminosity. The trend with luminosity is broadly con- 
sistent with the anti-correlation seen in local AGN, suggesting that 
this is the driving parameter rather than redshift. 

Unlike local AGN, however, we find evidence for an upturn 
in the variability at the highest luminosities (L x > 5 x 10 44 erg 
s~ 4 ). We speculate on possible explanations for this upturn. If lu- 
minosity is the driving parameter, it could reflect the turn over in 
the mass spectrum of active black holes. This could also be caused 
by redshift evolution, e.g. if luminous, high redshift QSOs contain 
smaller black holes or accrete more efficiently. 

Finally, we find evidence for X-ray variability in an object 
classified as a narrow emission-line galaxy, suggesting the presence 
of AGN activity. 
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In principle we can correct for each of these effects in order to 
compare all QSOs within the same frequency range. We therefore 
normalise all values of oq as follows: 



norm 



(1 + zf 2 ° x (3 x (Ti I week)~ x o d Q 



(A5) 



When combining an ensemble of QSOs, this correction can 
be included in Equation 7 (simply replace <tq with a'q ta as given 



APPENDIX A: CORRECTING FOR OBSERVATION 
LENGTH, BINNING AND TIME DILATION 

Detailed Fourier analyses of local AGN have revealed larger ampli- 
tude fluctuations at low frequencies, parameterised by power spec- 
tra of the form: 



p(f) oc r 



(Al) 



where a ~ 1.5. Thus the amplitude of variability will depend 
on the exact range of timescales that we sample. If we sample fre- 
quencies in the range /i < / < /2, the intrinsic variance will relate 
to the power spectrum as follows (a / 1): 



h 



P(f)df 



/2 



/l — o. rl — a 

1 - J 2 



(A2) 
(A3) 



One can use this to correct for different sampling by normal- 
ising all QSOs to the same frequency range. We adopt a low fre- 
quency bound /i corresponding to a timescale of 1 week (1.65 x 
10~ 6 Hz). This reflects the typical observation length in our QSOs. 
We arbitrarily choose /2 = 1CT 3 Hz, corresponding to the length 
of the shortest bins (~ 1000s). 

In order to normalise to these frequencies, 3 distinct correc- 
tions must be applied to the amplitudes (a) derived in Section 3. 
The largest of these arises because the 7 fields are of differing ob- 
servation lengths, varying from ~ 2 days for the QSF3 field to a 
maximum of ~ 14 days for the F864 field. We expect larger ampli- 
tude variations for longer integration times. For a total observation 
length T, assuming perfect sampling and /i <g /2, Equation A2 
leads to: 



(cr) OCT 



(l-a)/2 



(A4) 



A second bias arises because of the irregular gaps in the data 
and from the re-binning described in Section 2. In particular, the 
bins in many fainter QSOs are merged, which will reduce the ob- 
served variance by averaging out high frequency variability. Quan- 
tifying these effects would be exceptionally difficult analytically 
because of the irregular sampling. We therefore carry out simu- 
lations for each QSO in order to quantify this behaviour. This is 
described in detail in Appendix B. For the faintest QSOs with only 
2 widely separated bins, we find that these effects produce an aver- 
age reduction in <tq by a factor (3 = 1.34 compared to a perfectly 
sampled light curve. 

Finally a correction may also be applied for the effect of cos- 
mological time dilation. If we sample frequencies in a light curve 
in the range fi < f < fi then the intrinsic standard deviation in 
the light curve will be given by Equation A3. If instead we observe 
the redshifted frequencies /' = //(l + z) then the observed a will 
be reduced by approximately (1 + z)~ 0,25 . Thus the amplitude we 
have measured for QSOs of redshift z = 2 should be increased by 
a factor ~ 1.3 for comparison with QSOs at zero redshift. 



QSO). Since these three corrections rely on an assumed form for 
the power spectrum, which may not be entirely valid for the general 
QSO population, we also display results based on the raw, uncor- 
rected data for comparison. 



APPENDIX B: SIMULATING LIGHT CURVES 

A time process x(t) may be represented by its Fourier integral: 



x(t) 



X{f)e^ st df 



(Bl) 



The power spectrum of x(t) is defined as the modulus squared 
of the inverse Fourier transform (neglecting factors of 27r): 

2 



P(f) 



s 

J — t 



x(t)e 



-2nift 



dt 



\X{f)\ 2 



(B2) 



Rather than assuming random variability (ie. white noise) for 
the QSOs, we assume that all QSOs show the same power spectra 
as local Seyfert galaxies, ie. P(f) oc / _1 ' 6 (McHardy & Czerny 
1987, Lawrence & Papadakis 1993). A general theoretical light 
curve can then be constructed to satisfy this spectral form: 



x(t) 



f 

J — ( 



\P(f)\ 



1/2 e i[2nft- 



4>(f)] 



df 



(B3) 



where the real function <j>{f) represents the phase information 
lost by the power spectrum method. In a stochastic random pro- 
cess governed by a large number of variables we expect no relation 
between the phase and frequency of Fourier components, and in- 
deed this has been shown for a number of EXOSAT AGN (Krolik, 
Done & Madejski 1993). The light curves were therefore simulated 
using a random phase </>(/). Sampling x(t) over discrete time in- 
tervals tk, the Fourier integral collapses into the discrete Fourier 
transform and we obtain: 



x(tk) 



|P(/)| 1/2 COs[27T/t fc -<K/)] 



(B4) 



The lowest contributing frequencies are chosen at 10 Hz. 
The high frequency cut off (the Nyquist frequency) was chosen at 
an arbitrarily high value of 10~ 2 Hz, close to the limit of power 
spectrum observations (Nandra et al 1997). To obtain the required 
amplitude of variability, the light curves are scaled and a flat d.c. 
component added to give a total fractional variance of the required 
amplitude. The resulting light curve is then corrected to allow for 
the redshift of the QSO (viz. all times multiplied by 1 + z) and a 
random section is sampled and rebinned in exactly the same way as 
the real data. 1000 such simulations are carried out for each QSO 
in order to determine the correction factor, (3, which gives the mean 
change in the intrinsic standard deviation a as a result of the bin- 
ning structure: 



(3 = <Ti, 



cf (&binned) 



(B5) 
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For the 86 QSOs in our sample, values of this correction factor 
lie in the range 1 < f3 < 1.34 with the largest values for the faintest 
QSOs with only 2 widely spaced temporal bins. 



APPENDIX C: THE EVOLUTION OF BLACK HOLE 
MASS IN A LONG-LIVED QSO MODEL 

Averaged over an accretion timescale, the luminosity of a QSO re- 
lates to its accretion rate via: 



L(z) = eMc 2 



5.7X10M W 



M e yr- 



(Cl) 



where e is the accretion efficiency, taken to be ~ 10% after 
Rees (1984), and L is the bolometric luminosity emitted across the 
entire electromagnetic spectrum. For a typical radio quiet quasar 
the bolometric correction from 0.1 — 2 keV is approximately a fac- 
tor of 12 from gamma-ray to far infra-red wavelengths. 

In the long-lived QSO model, the evolution of the QSO lumi- 
nosity function (e.g. Boyle et al 1994) directly describes the evolu- 
tion of individual QSOs: 



Loc (1 + 2)" 



z < 2 



(C2) 



At higher redshifts, studies of optical (and radio) quasars sug- 
gest that the luminosity function is roughly flat during the epoch 
2 < z < 3, with a rapid decline beyond redshift z = 3 (Schmidt, 
Schneider & Gunn 1995). For the purposes of a long-lived model, 
we therefore assume QSOs began accreting at z = 3 and estimate 
the initial mass of the black hole by assuming that the system began 
by radiating at the Eddington luminosity (Rees 1984). 



1 o„ in38 M(g = 3) _x 

1.3 x 10 — — -ergs 



M: 



o 



(C3) 



In order to relate the evolution in time with the evolution in 
redshift, we use the standard relation from a Robertson-Walker 
metric (A = 0, q a = 1/2): 

t = 2/3i2'~ 1 (l + z)~ 3/2 



(C4) 



From the above, one can obtain the evolution of mass with 
redshift: 
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Figure CI. Showing the predicted growth in black hole mass from z = 
until the present day, assuming a long-lived QSO evolution model, with 
initial mass given by the Eddington limit. 



— — — (C5) 
dz dt ' dz 

= =^m(l + z)- 5/2 (C6) 
ecr 
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27ec 2 

We can integrate this relation to obtain the mass as a function 
of z, where we fix the constant of integration using Equation C3. 
This is shown in Figure CI. For the z < 2 epoch, this evolution is 
parameterised as follows: 

M (*) = ipr [l3 + 0.7£- 1 fc B - 1 (5.2-(l + *) 3 / 2 )] (C8) 
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